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Abstract. The three-state Potts field theory in two dimensions with thermal and 
magnetic perturbations provides the simplest model of confinement allowing for both 
mesons and baryons, as well as for an extended phase with deconfincd quarks. We study 
numerically the evolution of the mass spectrum of this model over its whole parameter 
range, obtaining a pattern of confinement, particle decay and phase transitions which 
confirms recent predictions. 
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1. Introduction 

The g-state Potts field theory describes the universality class associated with the 
spontaneous breaking of the permutation symmetry of q colors. In two spacetime 
dimensions the theory is defined for q < 4. The transition between the disordered 
(symmetric) phase and the ordered (broken symmetry) phase occurs as the temperature is 
varied. If an external magnetic field is allowed to be switched on, then the theory provides 
a model of confinement with several interesting features pQ. Below critical temperature 
and in absence of the field, the particle spectrum contains kinks, which interpolate between 
the q degenerate vacua. The magnetic field breaks (at least partially) the degeneracy of the 
vacua and those kinks which start or end on vacua that become false get removed from the 
spectrum. Those multikink configurations which consist of such kinks but start and end 
on vacua that remain true can give rise to new particles, which will be bound states of the 
confined kinks. The role of the kinks is analogous to that of quarks in chromodynamics, 
and the kink bound states arising play a role analogous to that of mesons and baryons, 
therefore the kinks and the bound states are sometimes referred to as quarks, mesons and 
baryons. 

The mechanism of confinement through the breaking of the degeneracy of discrete 
vacua is quite general in two dimensions. The Ising case (corresponding to q — 2) provides 
the simplest example and was first studied in [2] (see [3] for more references). In this as 
in other cases [H El E] , the single-component order parameter leads to bound states made 
up of a kink and an anti-kink, i.e. only mesonic particles are present. The q = 3 Potts 
model provides the simplest example allowing also for baryonic particles (made up of three 
kinks), as well as for an extended phase in the parameter space in which some kinks are 
deconfined. Moreover, in this case the renormalization group trajectories originate from 
a non-trivial fixed point. 

A qualitative characterization of the evolution of the mass spectrum in the three- 
and four-state Potts models for generic values of the temperature and of a magnetic field 
chosen to act on a single color was given in [Tj. Our aim in the present paper is to verify 
that picture in the q = 3 case by numerical calculations. We use the method called 
the truncated conformal space approach (TCSA) [7j, which is suitable for studying mass 
spectra. Its application requires that the theory is formulated as a perturbation of a 
conformal field theory. In the q = 3 case this is the -D4 minimal conformal model, and the 
fact that this is a 'non-diagonal' minimal model makes the situation technically somewhat 
more complicated compared to the usual applications of the TCSA. Interestingly, our 
investigation of the off-critical region allows us to fix also some previously unknown data 
concerning the conformal point. 

The paper is organized as follows. In section [2] we introduce the model and discuss 
the implementation of the TCSA, before presenting the results of the numerical study 
in section [3J Section 0] is devoted to the comparison with analytic results that can be 
obtained in the case of weak magnetic field, while section [5] contains a few final remarks. 
The results concerning the conformal point are given in the appendix. 
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2. Model and numerical method 

The three-state Potts model on the lattice [H [9] is a generalization of the Ising model in 
which each site variable s(x) can take three different values (colors) 1, 2, 3. At temperature 
T and in the presence of an external magnetic field H acting only on the sites with one 
specific color, say 3, the reduced Hamiltonian reads as 



where the first sum is over nearest neighbours. At zero magnetic field this Hamiltonian 
is invariant under the group S3 of permutations of the three colors. If the magnetic field 
is nonzero, then TC is invariant only under the group S2 = Z2 of permutations of the first 
two colors. 

The system described by the model has two phases, an ordered one and a disordered 
one. While the disordered phase is characterized by a unique ground state which is 
invariant under the Z2 symmetry, in the ordered phase the system has two possible ground 
states, which are interchanged by the Z2 transformation. The phase boundary separating 
the two phases in the (T, H) plane has a first-order and a second-order part. The second- 
order part runs through the half-plane with negative magnetic field and relates the q = 3 
Potts critical point at H = to the q = 2 (Ising) critical point at H = —00, where 
the third color is unaccessible. The first-order transition takes place below the critical 
temperature T c at zero magnetic field, where three ground states, which are permuted by 
the S3 symmetry, coexist; in the first-order transition the system also exhibits spontaneous 
magnetization (the variables cr a (x) = 5 s ( x ^ a — |, a = 1,2,3 have non-zero expectation 



The field theoretical description of the scaling limit of the model ([T]) begins with 
the observation that the S'3-invariant critical point at (T,H) = (T c ,0), where the first- 
and second-order transitions meet, belongs to the universality class described by the -D4 
minimal model of two-dimensional conformal field theory [10j QTl [12] . The corresponding 
value c = 4/5 of the central charge is the lowest one for which two different modular 
invariant partition functions can be obtained on a torus with periodic boundary condition 
|12j . The first realization, known as A 5 , contains, with multiplicity one, all the primary 
operators <fi whose conformal weight A^ appears in the Kac table, table [lj and describes 
the tetracritical point of the Ising model. The second realization, known as -D4, contains 
only a subset of the primaries in table (TJ some of them with multiplicity 2, as expected 
for the 3-state Potts model, which has a two-component order parameter (X) a =i o~ a = 0). 
The primaries entering the description of the 3-state Potts model are listed in table [2] 
together with their transformation properties under the group S3, which is the semidirect 
product of Z 2 and Z 3 ; contrary to the situation for the A 5 case, some operators in tabled 
have a non-zero spin s = A — A. 

The 3-state Potts field theory describing the scaling limit of (CQ) is obtained by 
perturbing the -D4 conformal theory with the leading thermal operator e and the leading 
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Table 1. Kac table for c = 4/5. 
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Table 2. Primary fields of D4 and their transformation properties. 
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Z2-even magnetic operator a + = (ai + er 2 )/\/2. This leads to the Euclidean action 

j #*<*)->> J P) 

where Ad 4 is the conformal part, r ~ mass 2_2Ae is related to the deviation of the 
temperature from its critical value and h ~ mass 2_2A<T + is proportional to the magnetic 
field. The field theory (j2J) describes a family of renormalization group trajectories 
flowing out of the origin in the (r, h) plane. They are conveniently parameterized by 
the dimensionless combinations 

T T 
V± = ( ±/l )(2-2A £ )/(2-2A CT+ ) = ( ± ^)9/14 ' ( 3 ) 

where the positive (negative) sign applies when h is positive (negative). 

As a (l+l)-dimensional theory, describes relativistic particles whose mass 
spectrum, measured in units of the lightest mass, is a function of rj± only. In order 
to study numerically the evolution of this spectrum we resort to the method known as 
truncated conformal space approach (TCSA) [7] . This involves first of all considering the 
theory on a cylinder of circumference R, in such a way that the Hamiltonian operator 
takes the form 

H =J[( L o + L o- ^) +rj^ e(x,0)dx-h^ a + (x,0)dx, (4) 

where the conformal part is expressed in terms of the zero index generators Lq and Lq 
of the chiral Virasoro algebras with central charge c = 4/5. At the conformal point the 
eigenvalues of Lq + Lq are simply the scaling dimensions A + A of the operators, so the 
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conformal space of states coincides with the operator space. Also the matrix elements 
of the perturbing operators e and cr + on this space of states can be computed exactly, 
since they are related to the structure constants of the conformal operator product 
expansion as 



In this way the eigenvalues of (j4j) can be determined by numerical diagonalization of the 
Hamiltonian matrix on a finite dimensional subspace of the conformal space. Obviously, 
the truncation of the space of states puts upper bounds on the number of levels and on 
the values of R which are numerically accessible. Nevertheless, with a number of states 
which is very reasonable from the point of view of computer time, it is normally possible 
to obtain a sufficiently accurate description of many energy levels in the lower part of the 
spectrum, up to values of R for which the asymptotic regime relevant for the theory on 
the plane is already visible. 

After its appearance [7j, the TCSA has been widely used, including for the study of 
multiple perturbations [13J, non-minimal models [14J and theories with boundaries [15] 
(see e.g. [16] for additional references). Its use for the case at hand would not present 
special difficulties were not for the following point. Although minimal (i.e. containing a 
finite number of operator families), the conformal theory we start from does not belong 
to the series of minimal models (known as A-series) for which the structure constants 
appearing in are known completely from the work of Dotsenko and Fateev [17] . The 
three-state Potts conformal point, as already said, belongs to the D-series of minimal 
models, and for this series not all the structure constants are completely known; more 
precisely, some of them are known only up to signs [18j US]- Of course, these signs are 
themselves essential for the determination of the energy spectra through the TCSA. As a 
matter of fact, we managed to determine them (at least those needed for our purposes) 
requiring that the output of the TCSA is physically meaningful. The results are given in 
the appendix. 

3. Evolution of the particle spectrum 

We used the TCSA to follow the evolution of the lower part of the spectrum of the 
Hamiltonian operator (j3j) in the whole range of the parameters r]±. More precisely, we 
numerically determined the spectrum at 80 points along the curve r 14 + \h\ 9 = constant 
shown in figure [TJ each renormalization group trajectory flowing out of the origin in the 
(r, h) plane and corresponding to a specific value of r/ + or r?_ intercepts such a closed 
curve at a specific point, in such a way that making a round trip along the curve amounts 
to spanning all the trajectories. Since mass ratios are fixed along a given trajectory, it is 
sufficient to determine them at a single point. 

For the study of the mass spectrum it is sufficient to consider the subspace of states 
with vanishing total momentum on the cylinder, namely the operators with zero spirQ. 

1 Notice that also the primaries with non-zero spin in table [5] contribute scalar descendants. 
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Figure 1. Phase diagram of the field theory |Q|). The thick lines correspond to the 
first-order (continuous line) and second-order (dashed line) phase transitions. Dots on 
the closed path r 14 + |/i| 9 = constant mark the points corresponding to the spectra shown 
in figures \^p,-h. The ordered phase corresponds to the dotted part of the path. 



Moreover, one can treat the parts of the Hamiltonian which are even and odd under the 
Z2 symmetry separately. We included in our numerical calculations conformal states up 
to level 8 of descendance within the conformal family associated to each primary operator. 
This corresponds to 2426 states in the even sector and 1829 states in the odd sector. 

Examples of the results we obtained for the spectrum at different values of rj± are 
given in figures [2^i-h. They show, as a function of the cylinder circumference R, the 
energy differences Ei — Eq, where Ei is the energy of the i-th level, with the ground state 
corresponding to i — 0. As a consequence, in these figures the ground state coincides with 
the horizontal axis. Due to the absence of phase transitions on the cylinder, the ground 
state is always unique at finite R. When other levels are seen to approach the horizontal 
axis as R increases, the spectrum refers to values of r and h corresponding to the ordered 
phase (with degenerate ground states) on the plane (R = 00). 

The energy levels behave as 1/R in the conformal limit R — > (see (jl])) and approach 
constant values as R — > 00. These asymptotic constants determine the particle spectrum 
of the theory on the plane, the one we are interested in in this paper. The lowest nonzero 
asymptotic value determines the mass of the lightest particle (or multi-particle) state 
with quantum numbers compatible with the boundary conditions chosen on the cylinder 
geometry. In this paper we always work with periodic boundary conditions, which select 
states with zero topological charge, and determine the mass of the kinks (which are 
topologically charged) in the ordered phase from that of the kink-antikink states. 

If mi is the mass of the lightest particle, the continuous part of the spectrum on the 
plane starts at 2m!. On the cylinder, the continuum breaks into infinitely many discrete 
levels which become dense when R — > 00: the lowest one (the threshold) corresponds to 
a pair of particles at rest, the others to a pair with momenta p and —p. Of course, in 
our truncated numerical approach, only a finite number of these 'momentum lines' are 
visible. In general, the theory on the plane possesses several stable particles with masses 
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(e) 7? ± = +oo (f) T] + = 1 

mj, and the total spectrum is made up of all the states resulting from the combination of 
these particles. 

With these remarks in mind, it is not difficult to see that the numerical results indeed 
confirm the evolution of the particle spectrum predicted in pQ. Referring the reader to 
that paper for the detailed arguments at the origin of the prediction, we recall here the 
main features and how they appear in the TCSA. 
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(g) V+ = -0.1 (h) V+ = -1.07 

Figure 2. The low lying energy differences Ei — Eq as functions of R at the points a-h 
along the closed path of figure [7J Even levels are in red; odd levels are in green with 
dashed lines. Units refer to a fixed mass scale. 



3.1. The ordered phase 

We start our discussion at low temperature (r < 0) and zero magnetic field, where the 
model is on the first-order phase transition line. Here the system on the plane exhibits 
a spontaneous breaking of the S 3 symmetry and has three degenerate ground states |0 a ), 
a = 1, 2, 3. The spectrum of elementary excitations is known [20] to contain only the kinks 
K a/ 3, a 7^ P, of equal mass m, which interpolate between the ground states a and (3. The 
finite volume spectrum calculated by the TCSA is shown in figure [2£i. Two superimposed 
levels are clearly seen to approach the horizontal axis to produce the asymptotic triple 
ground state degeneracy. The splitting of the vacua is due to tunneling effects and, 
up to truncation errors, decreases exponentially at large R. The unique finite volume 
ground state belongs to the Z 2 -even sector; it becomes |0 3 ) asymptotically, while the 
other two levels go into |0i) ± | O2) - The single-kink states do not appear due to the 
periodic boundary conditions, but two-kink states K a pKp a as well as three-kink states 
KapKp^Krya are clearly visible. 

When, starting from this situation, a magnetic field of the form specified in (J2J) is 
switched on, the vacuum | O3) is no longer degenerate with the other two: in particular, 
for h < its energy is higher and it becomes a false vacuum. In the finite volume this 
translates into an energy level linearly increasing with R, a feature perfectly visible in 
figure [2b. Since finite energy excitations cannot start from or end on a false vacuum, the 
kinks Kis, K23 and their antikinks get removed from the spectrum, while Kyi and K21 
survive. The multi-kink states in which 7 = 3 enters only as an internal index preserve 
a finite energy in the negative field. This is the case, in particular, of the state KisK 32 , 
which however cannot remain a two- kink excitation at h < 0, since the false nature of 
the intermediate vacuum prevents the two kinks from moving arbitrarily far apart. They 
are instead confined into the topologically charged 'mesons' irfy , with n a positive integer 
indexing mesons of increasing mass. While the mesonic spectrum becomes dense above 2m 
as h — > 0~, the number of stable mesons is expected to decrease as one moves away from 
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Table 3. Multiplicities of low lying states of zero spin in the representations of the 
conformal algebra entering the Ising model; n is displayed in the top row. 

the first-order transition line. Indeed, since the mass splittings increase in the process, 
more and more mesons become unstable crossing the decay threshold 3m, above which 
the decay channel ~* K12K21K12 opens. 

The pattern of the levels visible in figure [2b is consistent with these expectations. 
Because of the periodic boundary conditions, the topologically charged particles K\i and 
7T^ ) and their anti-particles appear only through neutral combinations!. So the mass 
gap corresponds to the KyiK^x threshold and the levels converging at 3m to the tt^-^i 
threshold. Figure [2b shows that at i]_ = the theory still is in the ordered phase 
(the single-particle excitations do not appear, so they are kinks). The absence of lines 
converging at three times the mass of the kinks now indicates that no stable topologically 
charged meson is left so far away from the first-order transition. The non-monotonic 
behaviour of the lowest lines in the even sector is a residual manifestation of the false 
vacuum visible in figure [2b, which is now very unstable. 

3.2. The second- order phase transition 

The reduction of the mass gap from figure [2b to figure [2b is consistent with the expectation 
that the ordered phase ends at a critical value nf_ where the kink mass vanishes and a 
second-order transition of Ising type takes place. The value n c _ had been naively set to 
in pp. There is, however, no symmetry (like duality in the tricritical Ising model) which 
selects this value and it was observed in [22] that, since the Potts Curie temperature is 
larger at q = 2 than at q = 3, nf_ is more likely positive. This latter expectation is indeed 
confirmed by the TCSA, which neatly shows a phase transition at T]°_ — 0.14 and confirms 
that it corresponds to a crossover from q = 3 to q = 2 criticality. 

To see this point notice that the signature of a crossover to Ising criticality in the 
finite volume spectrum is that for sufficiently large values of R the energy levels behave 
again as 1/R, or more precisely as 

Ei-E ~-£(Ai + Ai), (6) 

where Aj + Aj are scaling dimensions of Ising operators. Since we are considering the 
zero momentum sector, and on the cylinder the momentum coincides with the spin, only 
those states with A = A have to be taken into account. The Ising model with periodic 
boundary conditions contains three representations of the conformal algebra with highest 

2 The topologically charged mesons, whose stability is specific of the Potts model with q = 3, were not 
mentioned in [T]. They have been recently discussed in [2"T] , 
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(a) even sector 
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Figure 3. ^{E, L — Eq) as functions of R for the lowest values of i in the (a) even and 
(b) odd sectors at the second- order phase transition point. Degeneracies predicted by the 
critical Ising model (Table\B$ are shown on the right hand side. 
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Figure 4. The lowest energy difference -^-(Ei — Eo) in the two sectors as a function of R 
at the second- order phase transition point. Data are shown with dots, the constants 2 A 
equal 1/8 and 1 predicted by the critical Ising model are shown with continuous lines. 



weights A = 0, 1/2, 1/16. The scaling dimensions entering ([6]) are then 2(A + n), where 
n = 0,1,2,... labels the conformal levels. The multiplicities of the states with these 
dimensions are shown in table [3j 

Figure [3] shows ^(Ei — E ) as functions of R in the even and odd sectors at the 
point 77 _ = 0.14 (the same absolute ground state energy is subtracted in both sectors). 
According to ([6]) these functions should be approximately constant for sufficiently large 
values of R. The data are in good agreement with this expectation and also the location 
and the multiplicities of the lines agree well with the data of the critical Ising model. For 
the lowest energy difference in each sector the accuracy of the data can be appreciated in 
figure HI 

Similar studies of second-order phase transitions by means of the TCSA were done in 
[2"3| [21] . A second-order phase transition was also observed in the tricritical Ising model 
[5], however its type (which is known to be Ising) could not be established numerically 
because of large truncation errors. 
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3.3. The disordered phase 

When, starting from the first-order transition line, the magnetic field is switched on with 
a positive sign, | O3) remains as the unique true vacuum of the theory. We thus enter 
the disordered phase in which all kinks are confined inside topologically neutral mesons 
originating from K^iKx^-iiK^Kz^., and baryons originating from ^31^12-^23 ±-^32-^21-^13, 
where we have taken linear combinations with definite Z 2 -parity. We keep the notation of 
[lj and denote ir^ and tt^ the even and odd mesons, respectively, and p± the baryons 
with parity ±1. In all cases the positive integer n indexes particles with increasing mass, 
and again the spectra of these particles become dense (above 2m for the mesons and above 
3m for the baryons) in the limit h — > 0; as we move away from the first-order transition 
line more and more particles should cross the two-meson decay thresholds and become 
unstable. 

These expectations are fully consistent with the numerical spectrum shown in 
figure Eh, which exhibits a large number of lines corresponding to single-meson states and 
a doubly degenerate linearly rising line corresponding to the two false vacua. A linearly 
rising line-like feature with the same slope can also be seen somewhat higher, which is 
the signature of some topologically neutral configuration that is unstable because it is 
supported by the false ground stateqj. A line corresponding to the baryon p^ is also 
present, although it is not very clearly seen because of truncation errors and the presence 
of many other lines in its vicinity. 

At rj + = +00 the theory is integrable and is known [26], EJJ to possess a doublet of 
massive particles as the only single-particle excitations. By dualityB their mass coincides 
with the mass of the kinks at rj± = — 00. It is expected that the two lightest mesons 7Tq 
and 7Ti are the only products of kink confinement which do not decay as rj + increases, 
and that they become the above mentioned doublet at r] + = +00. This prediction is 
confirmed by the numerical calculations. The reduction in the number of stable single- 
particle states as r] + increases can be appreciated in figures [2g and [2f. In the latter only 
the two lightest mesons are left below threshold; their mass difference vanish at r] + = +00 
(figure [2b) and changes sign with the magnetic field (figure EH). 

3.4- Complete evolution 

The systematic collection and analysis of the TCSA spectra at regular small intervals 
along the closed path of figure [1] allows us to follow the evolution of the lightest particles 
of the theory over the whole range of the parameters r]±. The result is shown in figure 
and shows a remarkable agreement with the qualitative prediction given in figure 8 of pQ . 
The data show how the confinement of kinks at h > produces both mesons and baryons, 
and how all the mesons, excluding the two lightest, and all the baryons reach the decay 
thresholds as rj + increases. The mesons ir^ and ir^ are the only stable particles of the 
theory in a sufficiently wide range around T]± = +00, the point of enhanced (S3) symmetry 

3 See [25] for a discussion of (possibly stable) particles above threshold. 

4 See [1] for the relation between the scattering theories at r > and r < in zero field. 
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Figure 5. Masses of the lightest mesons Kq 1 ^, tt^, ^\ ^\ the lightest baryon 

p+ > the elementary kinks Kij and the lightest mesonic kinks ttQ , along the closed path 
of figure [JJ Even particles are shown in red, odd particles in green and the kinks in 
blue. Dashed lines show the stability thresholds for even particles (twice the mass of 
the lightest particle), odd particles (mass of the lightest even particle plus that of the 
lightest odd particle) and mesonic kinks (mass of plus that of ) in red, green and 
blue, respectively. The points a-h, having a correspondence in figures [JJ and\^ are also 
marked. 



where their mass trajectories cross. At h < the even meson 7Tq is no longer the lightest 
particle of the theory, and decays before the second-order transition point, where the mass 
of 7Ti vanishes, is reached. In the ordered phase the spectrum of stable particles contains 
the deconfmed elementary kinks K\2, and, close enough to the first-order transition 
point r]± = — oo, the mesonic kinks , ir^ . 

Of course, more and more particles emerge from the thresholds and become stable 
as r)± — > — oo. Figure includes only the five lightest topologically neutral mesons, the 
lightest baryon and the lightest mesonic kink. The behaviour of the curves at h = 
follows from the discussion in the next section. 



4. Weak magnetic field 

A number of analytic results can be obtained at h = 0, where the theory is integrable, 
and for small h, using perturbation theory around the integrable case. In this section we 
compare some of these results with the numerical ones that we obtain using the TCSA. 

At h = and, by duality, for both signs of r, the vacuum energy density and the 
mass of the elementary excitations can be written as 

e = ar 10/6 , m = br 5/6 , (7) 
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with a and b dimensionless constants whose exact values are known to b^fl [28] 

a = -5.856.., 6 = 4.504... (8) 

The measure of the slope of E and of the mass gap at large R in the TCSA gives a ~ —5 
and b ~ 4.5. We are able to determine a only with a low accuracy, a circumstance that is 
not unexpected: as a general experience absolute energies are given less precisely by the 
TCSA than mass gaps (energy differences). 

Consider now the high temperature phase, r > 0. The first-order correction to 
eo is proportional to the spontaneous magnetization, which is zero. The second-order 
correction reads 

Se = -L J d 2 x(a + (x)a + (0)> = dh 2 ^ 9 . (9) 

The susceptibility was computed in [29] (see also [3D], appendix B) within the two-particle 
approximation in the form factor approach, with the result d ~ —0.135, which compares 
very well with the result d ~ —0.133 that we obtain using the TCSA. 
Let us denote as 

5m ( p = ah + fih 2 + .. (10) 
the correction to the mass of the meson -k\ , i = 0,1. At first order one has [T31 [1] 

c i = -^<7r i (1) (0)| O - + (0)|7r«(0)> . (11) 

These matrix elements are known (see [30] and references therein) in a basis A, A in which 
vrf } = (A + (-1)VT)/V2, and read 

e = - ei = -0.7036.. r - 13/18 . (12) 

Our TCSA determination eo = —e.\ — — 0.71.. r~ 13 / 18 is in reasonable agreement with 
this exact result. As regards the second-order mass corrections, we obtain the numerical 
results /o ~ 0.09 r~ 41 / 18 and fi ~ 0.12 r~ 41 / 18 . These numbers probably have considerable 
errors, nevertheless it will be interesting to compare them with analytic results when these 
are available. 

Turning to the low temperature phase, the confinement of kinks at leading order in 
the magnetic field can be described within the non-relativistic framework pQ. The detailed 
analysis has been recently performed in [2T] for the mesonic spectrum, both at low energy 
and in the semiclassical limit, as a generalization of the Ising case. As regards the low 
energy part of the spectrum at h > 0, the linear confining potential in one dimension 
produces mesons with masses mf 1 which deviate from twice the kink mass by a term 
c ("0^2/3. constants are proportional to the zeros of the Airy function for the odd 
mesons (i = 1), and to the zeros of the derivative of the Airy function for the even mesons 

(i = o) m- 

A precise quantitative investigation of this region of the phase diagram in which the 
mesonic spectrum tends to a continuum requires a numerical accuracy higher than that of 

5 Throughout the paper we use the normalization of the operators in which (cj)(x)cj)(0)) — > |x| _4A * as 
\x\ ->0. 
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our TCSA data. Nevertheless, we found that for moderately large values of the magnetic 
field (rj + < — 1) the data for the masses of the five lightest mesons are described relatively 
well by the h 2 ^ 3 power law; the exponent that we could extract by fitting a power function 
to the data is approximately 0.7. As regards the prefactors, we obtain very approximately 
Ci/($P ~ 2, to be compared with the value close to 2.3 corresponding to the analytic 
prediction. 

5. Conclusion 

We studied the scaling limit of the two-dimensional three-state Potts model as a function 
of temperature and magnetic field. This has been done directly in the continuum limit, 
considering the perturbation of the D4 minimal model of conformal field theory by its 
leading thermal and magnetic operators. The magnetic perturbation was chosen to be 
coupled to one color only, in such a way to leave a residual permutation symmetry in the 
first two colors. We studied the particle spectrum of the theory in the whole plane of 
the coupling constants using the numerical method called the truncated conformal space 
approach, with the aim of verifying the qualitative predictions made in pQ . 

Our results confirm these predictions, clearly showing that kink confinement produces 
both mesonic and baryonic particles and that deconfined kinks survive within an extended 
ordered phase whose boundary contains a first-order and a second-order part. In 
particular, we determined numerically the location of the second-order transition, which 
is not fixed by symmetry. The evolution of the masses of the first few lightest mesonic 
and baryonic states has been followed through the whole parameter range of the theory, 
exhibiting the decay pattern and confirming that the two particles of the disordered phase 
in zero field are the two lightest mesons produced when breaking the ordered phase by a 
positive field. 

The implementation of the numerical method required the knowledge of the structure 
constants of the thermal and magnetic operators in the D4 conformal field theory. We 
determined the signs of the structure constants which were left undetermined in [TSJ [TH] . 
The integrability of the model at zero magnetic field makes it possible to obtain analytic 
results for the particle masses and for the vacuum energy density at zero and small 
magnetic fields. We compared our numerical results with such analytic results at high 
temperature, finding a good agreement. We also presented numerical results for the 
second-order corrections of the meson masses, which could be useful for comparison with 
future analytic results. We were able to partially confirm recent results of [21] for the 
meson masses at low temperature and small magnetic field. This region of the phase 
diagram is numerically more demanding because the mesonic spectrum becomes dense as 
the magnetic field approaches zero, so that a complete verification requires a numerical 
accuracy larger than that of the present study. Normally this can be achieved by increasing 
the truncation level. 
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Note added 



Following the appearance of this article as a preprint, G. Takacs informed us that complete 
results for the structure constants of the D-series minimal models were obtained in [3T] . 
We checked that, once the different conventions and normalizations are suitably taken 
into account, the structure constants of [31] for the _D 4 case exactly coincide with those 
that we give in the Appendix. 



Appendix 

The structure constants for the D-series unitary minimal models were studied in e.g. 
J. They can be written as 



Cyk — M ijk J D A . A . Ak D A . A . Ak , (13) 

where A», Aj are the left and right conformal weights of (pi, D A . AjAk and D AiA . Ak are 
structure constants of the corresponding A-series minimal model (for the D4 model this is 
the A 5 model), and are called relative structure constants. The constants D AiAjAk 
and D AjAjAk can be calculated using the formulae given in [T7]. For the _D 4 case the 
relative structure constants are real numbers on which the £3 symmetry and 

a\ = a 2 , Q\ = n 2 (14) 

impose some restrictions. The absolute values of these constants were found to be either 
or 1 or v2- The latter value occurs in the following cases: 

[ GTl (Tl j I -^^<T1<T2<T2 I 

\Ma 2(Tiai \ = \Mq 1(T2(T2 \ = iM^^nJ = \M air72 n 2 \ = y/2. (15) 

The signs of the relative structure constants were not completely fixed in [TS1 [T5] . 
It is useful to introduce the combinations 

which are even and odd, respectively, under Z2 transformations. The model has an even 
part consisting of the modules /, e, X, Y, cr+, Q + , and an odd part consisting of the 
modules o~_, Q_, J, J, W, W. The fields in the even part form a closed subalgebra, 
which is also contained by the A§ model. This means that the even part of the D4 model 
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coincides with a sector of the A§ model. The \/2 mentioned above is eliminated from the 
relative structure constants by using the even and odd combinations introduced above. 

Since in this paper we study the D 4 model with the e and a + perturbations, the 
relative structure constants that we actually need to know are M it j and M ia+ j. with 
j = e or o~+ is nonzero if both D^^^ and D^.^.^ k are nonzero, and the S3 symmetry 
and (fT4|) also allow a nonzero value. Then, the nonzero M ie j are 

M Iee , M eeI , Mxee, M etX , M XtY , M YtX , M n tn ..\l n .,n .-\lv.<n.. (18) 
for i and j in the even sector, and 

M a _ ea _ , M a _ en _ , M n _ ea _ , M j€ j, Mj tJ , M JeW , M WeJ , M JeWl M WeJ , (19) 
for i and j in the odd sector. The nonzero M ia+ j are 

My CT+0 . . M a+a+ i, \ /„.„.,. „ , . . f , . ,y. 

^si +( j + i, M^^y, M a+a+a+ , Mn +a+u+ , M a+a+ n + , (20) 
for i and j in the even sector, and 

Mj a+a _ , M Ja+n _ , M J(T+CT _ , M jCT+n _ , 

Mn^j, M a _ a+W , M a _ a+W , M a _ a+n _ , M n _ a+a _ , (21) 

for i and j in the odd sector. 

The signs of the constants listed in (fT8"|) . ffTQl . f[2"Ul) and fT2TT) are partially determined 
by the 5*3 symmetry and by (Tbil) . The constants (|T8l) and (l20i) are equal to 1 because the 
even sector coincides with a sector of the A 5 model. 

In order to find the signs of the constants f|T9|) and (121]) for the odd sector we used 
the empirical criterion that physically meaningful spectra should be obtained if the -D4 
conformal theory is perturbed by the e and a + fields. The TCSA calculations allow one 
to fix the signs of all constants (fT9|) and (1211) up to a trivial freedom that corresponds to 
sign changes of the basis vectors. Within one of these equivalent choices the result that 
we obtained is that the constants (Tl9l) are equal to 1; the constants (121]) are equal to — 1, 
with the following exceptions: 

M n _ a+J = M n _ a+J = M Ja+n _ = M Ja+n _ = 1. (22) 
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